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ABSTRACT 

The All Sky Automated Survey (ASAS) monitors bright stars (8 mag < V < 
14 mag) south of declination +28°. The ASAS Catalogue of Variable Stars (ACVS) 
presently contains 50,099 objects; among them are 2212 objects classified as RR Lyrae 
pulsating variables. We use ASAS photometric V band data to search for multiperi- 
odicity in those stars. We find that 73 of 1435 RRab stars and 49 of 756 RRc stars 
exhibit the Blazhko effect. We observe a deficiency of RRab Blazhko variables with 
main pulsation periods greater than 0.65 days. The Blazhko periods of RRc stars ex- 
hibit a strongly bimodal distribution. During our study we discovered the Blazhko 
effect with multiple periods in object ASAS 050747-3351.9 = SU Col. Blazhko peri- 
ods of 89.3 d and 65.8 d and a candidate of 29.5 d were identified with periodogram 
peaks near the first three harmonics of the main pulsation. These observations may 
inspire new models of the Blazhko effect, which has eluded a consistent theory since 
its discovery about one hundred years ago. Long term light curve changes were found 
in 29 stars. We also found 19 Galactic double mode pulsators (RRd), of which 4 are 
new discoveries, raising the number of ASAS discoveries of such objects to 16, out of 
27 known in the field of our Galaxy. 
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1 INTRODUCTION 

The discovery of RR Lyrae stars in star clusters (Bailey 
1902) led to some of the first studies of radial pulsations 
in stars (Shapley 1914, Eddington 1917) as well as distance 
determinations based on stellar variability (Shapley 1918). 
Recent surveys have searched for substructure in the Milky 
Way, attributed to the incorporation of satellite galaxies, 
via overdensities of RR Lyrae stars (Duffau et al. 2006). 
Their pulsation means they lie in the instability strip of the 
Hertzprung-Russell diagram, the extent (and metallicity de- 
pendence) of which may be evaluated based on a statistical 
sample of such stars. Finally, RR Lyrae stars are found to 
pulsate not just with a single period, but in multiple ra- 
dial and non-radial modes of different period, which gives 
constraints both to pulsation theory and to internal struc- 
ture calculations. The current study presents the photomet- 
ric signal of these multiperiodic phenomena in a large sample 
of RR Lyrae stars. 

In Section 2 we shortly describe ASAS experiment and 
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its database, as well as the RR Lyrae sample derived from 
the ASAS Catalogue of Variable Stars. Sections 3 and 4 
contain the method and results of multiperiodicity search, 
together with the discussion. Finally in Section 5 we sum- 
marize the results. 



2 ASAS DATA 

2.1 The All Sky Automated Survey 

The ASAS project is located in Las Campanas Observatory 
in Chile and currently consists of two small, wide-field tele- 
scopes (200/2.8) with 2K x 2K CCD cameras with 15/x pixel 
size from Apogee, observing in standard V and I filters. 
The main goal of the ASAS project is to monitor the whole 
available bright (8 mag < V < 14 mag) sky and search for 
variability. It has been observing south of declination +28° 
(almost 75% of the whole sky) since 2000, and a smaller part 
since 1997, covering the available sky every 1-3 nights. The 
project is ongoing, but has already yielded a few hundred 
photometric measurements per star. 

Analysis of the V band data taken over the last 
6 years resulted in the ASAS Catalog of Variable Stars 
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(ACVS), which presently contains 50,099 variable stars. 
The / band data is currently being processed and is 
not yet included in the catalog. ACVS is available both 
as an online database and a downloadable text file. 
For details on ASAS equipment and catalogues see Po- 
jmahski (1997, 1998, 2000, 2002, 2003), Pojmahski and 
Maciejewski (2004, 2005), Pojmahski, Pilecki and Szczy- 
giel (2005). All the ASAS data that is public domain, as 
well as an online catalogue, are available on the WWW: 
http:/ /www.astrouw. edu.pl/~gp/asas/asas. html 
http:/ /archive. princeton.edu/~asas/ 

2.2 RR Lyrae Sample in ACVS 

ACVS contains a sample of 2212 RR Lyrae stars: 1455 of 
RRab type and 757 of RRc type. In variable stars with am- 
biguous classification, we used those for which RR Lyr type 
is the most probable. Their distribution across the sky in 
galactic coordinates is shown in Figure 1. 

The classification procedure of variable stars in ASAS 
was based mainly on Fourier parameters of the light curve, 
JHK colours from 2MASS, and the object's period as de- 
scribed in Pojmahski (2002). The full catalogue consists of 
five parts created between 2002 and 2005, so the classifi- 
cation of the light curves as RR Lyrae stars was based on 
measurements taken from different time baselines. But we 
used the longest available baseline (for most objects between 
HJD 2451800 to HJD 2453800) to refine their periods, which 
was especially useful for objects from the oldest catalogues. 
We also used all available data points of reasonable qual- 
ity (i.e. with flags A and B) in the search for multiperiodic 
behaviour. 

Because the ASAS classification was automatic and al- 
lows multiple type assignment we decided to roughly inspect 
the light curves visually to remove variable stars of other 
types, which were mistaken for RR Lyrae stars. That left 
us with 1435 RRab and 756 RRc stars. It is important to 
mention that the inspection was very superficial and ex- 
cluded only obvious misclassifications among RRab objects, 
e.g. MIRA type variables. The light curves of RRc type are 
very similar to EC (i.e. eclipsing contact binary) type, we 
simply trust the automatic classification, accepting a high 
rate of misidentification. 

The distribution of periods of our RR Lyrae sample is 
shown in Figure 2. Figure 3 presents two histograms also 
based on the sample, one of magnitudes in maximum light 
and the other one of the amplitudes of the light curves. 



3 BLAZHKO EFFECT 

The shape of the light curves of many RR Lyrae stars 
changes in a nearly periodic fashion, with time-scales of be- 
tween about five and a thousand or more days; this is called 
the Blazhko effect (Blazhko 1907, Smith 1995). Since this 
shape change is equivalent to periodic phase and amplitude 
changes of the harmonics which make up the light curve, 
in the frequency domain it is identified through sidebands 
(additional peaks) in the periodogram near the harmonics of 
the main pulsation. Physically these peaks may correspond 
to excited non-radial modes, but the physical nature of the 
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Figure 1. Distribution of RR Lyrae stars from ASAS Catalogue 
of Variable Stars (ACVS) in galactic coordinates. The dashed line 
is the equatorial plane. 
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Figure 2. The distribution of periods of the RR Lyrae sample. 



Blazhko effect is still unsettled a century after its discov- 
ery. For some stars, which we label BL1, there is a single 
significant frequency near the main mode. For other stars, 
which we label BL2, two peaks flank the main frequency 
with approximately equal frequency spacing (Mizerski 2003, 
Alcock et al. 2003). There are also some stars with two ad- 
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Figure 3. The distribution of V magnitudes in maximum light 
and V amplitudes of the RR Lyrae sample. 
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ditional peaks which do not form an equidistant triplet with 
the main pulsation and others with more than two peaks. 

3.1 Identifying extra frequencies 

The method we use for discovering Blazhko stars is a slightly 
modified version of the method applied by Collinge et al. 
(2006) to Optical Gravitational Lensing Experiment data. 
That paper gives a full introduction to the method which 
we discuss more concisely here. 

We first subtract a Fourier series of best-fitting har- 
monics of the main pulsation period from the data, which 
is called prewhitening. The number of harmonics to fit was 
chosen by the method of Kovacs (2005) and depends on an 
empirical signal to noise level. We next run the CLEAN algo- 
rithm (Roberts et al., 1987) on the prewhitened light curves, 
using 200 iterations and a gain of 0.5, on the frequency band 
[0-6 day -1 ] with a frequency resolution of 0.125/T, where 
T is the overall time baseline of the observations for the 
given light curve. The peaks are sorted by amplitude, and 
the noise level is indicated by the 15th largest peak (Aiz). 
By a Monte-Carlo method, we determine how much bigger 
than that noise level a signal must be to classify as a de- 
tection. The Blazhko effect has a period of longer than ~ 5 
days, so we define the "Blazhko range" as the frequencies 
within 0.2 day -1 of, but further than 1/T from, the main 
pulsation frequency. 

BL1 stars are those with a single peak with amplitude 
A > 2.0*4i5 in the Blazhko range of the CLEANed spec- 
trum. BL2 stars have two frequencies within the Blazhko 
range, on opposite sides of the main frequency, equidistant 
to it to within 3.0/T, and both frequencies must have am- 
plitudes A > 1. 15^li5. We follow Alcock et al. (2003) in 
defining a class of period change (PC) stars that have one 
or more peaks within 1 /T of the main frequency with am- 
plitude^) A > 1.3.4i5. The frequency criterion for PC stars 
means that even if the phenomenon is periodic, we have not 
yet observed a full cycle; we have only detected a chang- 
ing light curve shape. We note that a long-term amplitude 
change can also give an extra frequency component near the 
main frequency, so the PC class is imprecisely named. 

During its many iterations, CLEAN produces somewhat 
biased amplitudes, so after selecting the significant frequen- 
cies we refine the fit parameters with a non-linear solver (the 
Levenberg-Marquardt algorithm; Press et al. 1989); the re- 
sulting frequencies and amplitudes are reported in Tables 
2, 3, and 4. The best-fitting model is subtracted from the 
light curve and the amplitude of the highest remaining peak 
within 0.2 day -1 of the main frequency is also reported in 
Tables 2, 3, and 4. We do not consider this peak significant; 
it is reported to quantify the level of the noise for complete- 
ness studies. 

We determined the amplitude thresholds cited above 
by creating 3 false catalogues, each of them based on the 
real data sets for each star. In these catalogues the magni- 
tude residuals from the pre-whitening step were randomly 
swapped among the observation times. We search for peaks 
in the false catalogues with the CLEAN algorithm and set 
the amplitude thresholds to allow an acceptable false alarm 
rate such that about 2% of the stars in each of the categories 
BL1, BL2, and PC are likely to be spurious. 

This method allowed us to determine unambiguous 



Table 1. Numbers in different categories of RR Lyrae stars ex- 
hibiting Blazhko-typc behaviour. See section 3.2 for details 
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threshold values in the BL2 and PC groups, but it par- 
tially failed for the BL1 category. It turned out that false 
catalogues very often contain a frequency in Blazhko range, 
whose amplitude is above the chosen noise level. This is 
probably a result of high noise in ASAS data. And while 
the probability to find BL2s in false catalogues is low as it 
requires two evenly spaced frequencies, it grows rapidly for 
BLls, where only one frequency is needed to create a false 
alarm. 

So for BLls we used a more conservative method of 
threshold determination. Each prewhitened light curve was 
separated into even-numbered and odd-numbered points 
in the chronologically ordered data set, making two light 
curves. Then the search for Blazhko frequency was per- 
formed in all three light curves — the original and subdivided 
ones. A frequency must appear at least in two of the three 
light curves to qualify as a true detection. This search was 
performed for each star in the real and the 3 false catalogues, 
which allowed us to determine an unambiguous threshold 
value, because in random realization of the light curve such 
coincidence is less likely to appear. 

It turned out that in a real catalog 22 BL1 stars had a 
Blazhko frequency found in all 3 light curve realizations and 
30 BL1 stars had it in two, while in false catalogues it never 
occurred in all 3. Those stars which had a BL1 frequency in 
all 3 light curves are marked with a "*" in the last column 
of Table 2. 



3.2 Objects exhibiting Blazhko effect 

The number of unique stars in each of our categories is (BL1, 
BL2, PC) = (52, 70, 29); see Table 1 for details. 

We see that the number of Blazhko stars, namely 122 
objects, is very low; we found the Blazhko effect in 5.1% of 
RRab stars and 6.5% of RRc stars, which is much lower than 
the number determined by Mizerski (2003) for the Galactic 
Bulge (25% RRab and 10% RRc) and Soszynski et al. (2003) 
for the LMC (15% RRab and 6% RRc). The result for the 
SMC by Soszynski et al. (2002) of 10% RRab and 10% RRc 
is somewhat closer to our result, but the discrepancy is still 
visible. The photometric precision of the data is a few per- 
cent, which suggests there is serious incompleteness for the 
Blazhko effect at low amplitudes. Figure 4 (left panel) shows 
the distribution of mean V magnitudes for Blazhko stars, 
which is somewhat shifted towards the bright end relative 
to all RR Lyrae stars (Figure 3), suggesting that detection 
fails for faint stars, for which data scatter is higher. 

We also compared our results with previous BL searches 
in the sky surveys of the Galaxy. In the recent paper on a 
similar data quality Northern Sky Variability Survey (Wils, 
Lloyd and Bernhard 2006) only 34 out of 785 (4.3%) stars 
exhibit the Blazhko effect, although the time span of the 
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Figure 4. The distribution of V magnitudes of Blazhko stars (left), Blazhko periods (middle) and Blazhko V amplitudes (right) in the 
RR Lyrae sample. 
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Figure 5. Main pulsation period vs. Blazhko period. The Blazhko 
periods of RRc stars exhibit a strongly bimodal distribution. 



data was just a year. In a Blazhko search in ASAS data done 
by Wils and Sodor (2005), only 43 objects were identified 
(~ 4%), though the search did not include the latest part of 
ACVS with objects having declinations 8 > and inspected 
only RRab variables. 

While the result in this paper is an improvement com- 
pared to numbers given by other Galactic field investiga- 
tions, the percentage of Blazhko stars is significantly lower 
than for other systems. What is even more striking, is that 
we found slightly more Blazhko object among RRc variables 
than among RRab which is opposite to the usual findings. 
We think that the result for the RRc group is real (the per- 
centage of Blazhko is similar to the one in previous papers, 
that is 5-10%), but for some reason we are missing Blazhko 
objects among RRab stars. It is possible that the RRabs 
with strong Blazhko effect have such scattered light curves 
that they are classified as different type variables by the 
automatic ASAS classification. 

Figure 5 plots the Blazhko period Pbl versus main pul- 
sation period P. It is interesting that there is a significant 
difference in behaviour of RRab (filled symbols) and RRc 



(open symbols) groups. RRab stars have a relatively uniform 
distribution in the logarithm of the Blazhko period, whereas 
RRc stars occupy either a short Pbl region around 10 days, 
or a long Pb l region, starting around 300 days with a con- 
densation roughly around 1500 days. That is, we did not 
observe Blazhko changes with periods above ~ 20 days and 
below ~ 300 days in RRc pulsators. The unusual concentra- 
tion around Pbl 1500 days, which is often close to the time 
span of the data, is adding to the peak near Ioq{Pbl) = 3 in 
Figure 4 (middle panel). This concentration may be a false 
signal due to overall bias and those object might in fact be 
just exhibiting long period changes (PC). But there are ob- 
jects having Pbl values of 300-500 days, which is 4-5 times 
shorter than the data span, thus the gap cannot be easily 
explained as a pile-up of PCs misclassified as BLs. 

Jurcsik et al. (2005) found a correlation between the 
pulsation period and the Blazhko period, such that RR 
Lyrae with periods P < 0.4 d can have Pbl of order of days, 
while those with longer pulsation periods such as P > 0.6 d 
always exhibit modulation with Pbl > 20 d. Our Figure 5 
does not support this statement. They also suggested that 
there is a "continuous transition" between BL1 and BL2 
type modulations. While we observe differences in Blazhko 
behaviour between RRab and RRc groups, which can be 
due to physical differences in pulsation mechanisms in these 
groups, we note that BL1 (circles) and BL2 (triangles) oc- 
cur in roughly constant proportion throughout Figure 5. 
This observation may suggest that both BL1 and BL2 ef- 
fects have the same origin, it is just that in case of BL1 we 
do not observe the second equidistant peak present in BL2 
category. 

Figure 5 also reveals that the number of BL stars falls 
rapidly for pulsation periods P > 0.65 d (log(P) > —0.2). 
We compared the period distributions for RRab stars ex- 
hibiting the Blazhko effect with those that do not. Figure 
6 is a plot of the cumulative distribution function of the 
main pulsation period for Blazhko RRab stars and for RRab 
stars whose light curves are not modulated. The maximum 
difference is 0.241466 which occurs at P = 0.552945 d and 
has a false alarm probability of 4.8 x 10~ 4 according to the 
Kolmogorov-Smirnov test, as computed with the KSTWO 
function from Press et al. (1989). This observation could be 
an important constraint for models of the Blazhko effect, 
but it may also indicate misclassification of RRab stars that 
are modulated. That is, light curves of BL RRab are not as 
well-defined as their non-Blazhko counterparts and so have 
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Figure 6. The cumulative distribution function of the main pul- 
sation period for RRab Blazhko stars and those RRab that are 
not modulated. The maximum difference in the curves is marked 
with a bar and is highly significant. The lower incidence of the 
Blazhko effect for stars with long pulsation periods could give 
fundamental insight to the physical cause of the Blazhko effect, 
but we suspect that they exist but dropped out of our sample due 
to misclassification (see text). 



different Fourier parameters, on which the classification of 
Pojmahski (2002) was based. The risk of misclassification 
is especially strong for long period BL RRab, which may 
appear as 8 Cephei first overtone pulsators if their fourth 
harmonic is reduced relative to the second harmonic (see 
Figure 6 of Pojmahski, 2002). If this explanation is true, it 
would help account for the low percentage of RRab BL stars 
with respect to other surveys, since there are some Blazhko 
stars that were not found because we accepted the automatic 
classification. 

On the other hand, this observation is more likely a 
confirmation of a previously recognized effect, that the mean 
main pulsation period of Blazhko RRab stars is lower than 
that of the unmodulated RRab stars (Preston 1964, Smith 
1981, Nemec 1985, Gloria 1990). This effect has started to be 
considered in models of the Blazhko effect (Stothers 2006), 
and the current study is the most statistically significant 
detection of the effect to date. 

We also checked if there is any dependence of Blazhko 
amplitude on pulsation period. We do not see any significant 
trend; the distribution seems to be uniform. The right panel 
in Figure 4 shows the distribution of modulation amplitudes. 
Similarly, Blazhko amplitude versus modulation period (Fig- 
ure 7) does not present any relation. Both long and short 
Blazhko periods can have an amplitude as small as 0.02 mag 
and as big as 0.2 mag. 

The detailed results of our searches are listed in Table 
2 for BL1 stars, Table 3 for BL2 and Table 4 for PC objects. 

Table 2 contains magnitude, periods and amplitudes for 
BL1 stars. The main pulsation period is P = f^ 1 , where 
fo is the main pulsation frequency, and Blazhko period is 
defined as Pbl = (/i — /o) _1 , where fi is the frequency 
component caused by the Blazhko effect. Thus the negative 
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Figure 7. Blazhko amplitude versus Blazhko period. 

value of Pbl means that the additional frequency is smaller 
than the main pulsation frequency. This occurs in 10 out 
of 20 cases of RRc type and only 7 out of 32 of RRab, 
consistent with previous statistics (eg. Moskalik and Poretti, 
2003, Mizerski 2003) 

Tables 3 contains magnitude, periods and amplitudes 
for BL2 stars. The main pulsation period is also P = fg 1 , 
where fo is the main pulsation frequency. Now the Blazhko 
period is defined as Pbl = (/-,+ — /o) _1 where /_,+ means 
we choose between Blazhko frequencies (/_ and /+, with 
/_ < fo < /+) the one with higher amplitude. The negative 
value of Pbl means that the Blazhko peak of higher ampli- 
tude has a frequency lower than the main pulsation. This 
corresponds to the ratio of A-/A+ > 1. We see that this 
case occurs in 13 out of 29 objects of RRc type and 19 out 
of 41 of RRab. 

Finally, Table 4 contains magnitude, period, and am- 
plitude information for PC stars. The designations are the 
same as in Table 2. For RRab, equal numbers of stars have 
the modulation frequency lying on left and right side of the 
main pulsation frequency, which suggests that PC stars may 
be physically distinct from BL1, and not simply BLls with 
periods longer than the dataset. 

Two stars (one RRab and one RRc) have two pairs of 
equidistant frequencies in the prewhitened spectrum. Their 
properties are listed in Table 5. The RRab star, SU Col, 
is discussed in the next section. For the RRc star, LS Her, 
these pairs of BL2 peaks are close to one another, probably 
indicating that the Blazhko effect has an unstable period. 

3.3 A star with a rich Blazhko spectrum 

In the course of this analysis we identified one star (ASAS 
050747-3351.9 = SU Col) that has two clearly distinct pairs 
of BL2 peaks. In order to look for additional extra fre- 
quencies, we first folded the raw light curve and deter- 
mined, based on the scatter, that a 6 harmonic model would 
be appropriate. After removing those, we ran the CLEAN 
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Table 2. Objects exhibiting BL1 effect. Ao is the pcak-to-pcak amplitude of variation as provided by ACVS and A\ is the peak-to-pcak 
amplitude of Blazhko modulation. See section 3.2 for details. 



ASAS ID V P P BL A a A x Ax/Aq A noise Type Other ID 

[mag] [days] [days] [mag] [mag] 



001141-0144.9 


11.76 





.5297456 


154 


.53 





.97 





12 





12 





03 


RRAB 


R.Y Psc 




003338-1529.2 


11.09 





.5737309 


-255 


.50 





.64 





.11 


0. 


18 





04 


RRAB 


RX Cct 


* 


003706-4317.7 


13.29 





.6275343 


187 


.12 


1 


.31 





15 


0. 


11 





.04 


RRAB 


- 




011831-6755.1 


11.42 





.4057948 


1748 


86 





.49 





11 


0. 


23 





.04 


R.R.C 


AM Tuc 


* 


013140-4957.3 


12.12 





.4604329 


40 


.17 


1 


.10 





.18 


0. 


17 





02 


RRAB 


NSV00539 


* 


020752-2651.9 


9.48 





.4954335 


34 


.88 


1 


13 





.09 


0. 


08 





.02 


RRAB 


SS For 




022637-4119.7 


10.08 





.2941932 


357 


.94 





.13 





.03 


0. 


26 





.01 


RRC: 


- 


* 


025021-6415.7 


12.80 





.5724975 


-79 


.81 


1 


.29 





.16 


0. 


13 





04 


RRAB/DCEP-FO 


RV Hor 


* 


031113-2629.0 


11.43 
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31 


.81 


1 


.02 





12 


0. 


11 





02 


RRAB 


RX For 




033108+0713.4 


10.64 





.5281963 


-1442 


.79 





.18 





.06 


0. 


32 





.02 


RRAB/EC/ESD 


- 




053022-3234.8 


11.77 





.2331114 


364 


.60 





.20 





.05 


0. 


24 





.01 


RRC/EC 


- 




054230-1622.9 


11.85 





.5389135 


22 


.84 





.86 





.09 


0. 


10 





.03 


RRAB 


- 


* 


055322-5417.9 


12.82 





.2452638 


381 


.58 





.56 





.17 


0. 


31 





.05 


RRC/EC 


- 


* 


061401-6128.4 


9.31 





.4857445 


-117 


.90 





.46 





.02 


0. 


04 





.01 


RRAB 


ST Pic 


* 


062326+0005.8 


11.96 





.5501317 


37 


.24 





.63 





15 


0. 


23 





04 


RRAB 


- 




064615-4319.2 


13.05 





.3186021 


1639 


.61 





.71 





20 


0. 


28 





.07 


R.R.C 


- 




070854+1919.7 


8.69 





.7789192 


7 


.92 





.10 





.05 


0. 


45 





02 


RRAB 


- 


* 


071549-4405.3 


13.04 





.3144609 


-488 


.40 





.63 





.20 


0. 


32 





.05 


RRC 


- 




081933-2358.2 


10.45 





.2856671 


-8 


.10 





.29 





13 


0. 


47 





06 


RRC/EC 


- 


* 


085254-0300.3 


12.42 





.2669022 


-11 


.80 





51 





.23 


0. 


46 





.09 


RRC/EC/ESD 


- 


* 


090900-0410.4 


10.68 





.3032613 


-8 


.52 





.41 





.09 


0. 


22 





04 


R.R.C 


- 


* 


093731-1816.2 


13.05 





.5209175 


87 


.73 


1 


.23 





24 


0. 


20 





.06 


RRAB 


- 




101200+1921.9 


10.52 





.4826394 


-1141 


.03 





.85 





02 


0. 


03 





.13 


RRAB 


- 




123812-4422.5 


12.69 





.5235490 


1307 


70 


1 


.28 





22 


0. 


17 





.06 


RRAB: 


- 




124805-0820.8 


11.98 





.5105287 


65 


.69 





.88 





14 


0. 


16 





.04 


RRAB 


- 




135740-1202.3 


12.18 





.2671226 


1189 


91 





.44 





14 


0. 


32 





.03 


R.R.C 


- 




135813-4215.1 


12.43 





.5231816 


-146 


.01 





.84 





13 


0. 


Hi 





.04 


RRAB 


- 




141025-2244.8 


12.50 





.6398808 


1556 


66 





.82 





.11 


0. 


13 





03 


RRAB 


- 




144154-0324.7 


11.40 





.2293674 


-5 


.65 





.32 





.06 


0. 


18 





.01 


RRC/EC 


- 


* 


145315-1435.9 


12.43 





.5400738 


41 


.77 





.94 





.21 


0. 


23 





.05 


RRAB 


- 




155553-4041.7 


11.54 





.5819812 


48 


.82 





.72 





13 


0. 


18 





.05 


RRAB 


NSV07330 


* 


170223-2422.0 


11.34 





.4613693 


22 


.18 





.40 





.05 


0. 


14 





.01 


RRAB 


- 




172721-5305.9 


12.30 





.4354330 


58 


.66 


1 


.25 





15 


0. 


12 





.03 


RRAB 






180023-7026.5 


12.08 





.3556146 


1162 


.79 





.46 





.06 


0. 


13 





.02 


R.R.C 






181215-5206.9 


12.59 





.8375462 


-5 


.22 





.67 





.13 


0. 


19 





01 


RRAB 






185719-6321.4 


12.28 





.4120017 


61 


.39 





.96 





.12 


0. 


13 





.03 


RRAB 






194502+2434.2 


11.72 





.8458661 


37 


.56 





.30 





.06 


0. 


20 





.01 


RRAB/EC/ESD 






200431-5352.3 


10.95 





.3002402 


10 


.82 





.31 





08 


0. 


26 





.04 


RRC 






200556-0830.9 


12.75 





.4381966 


192 


.20 


1 


12 





18 


0. 


16 





.04 


RRAB 


KM Aql 




202746-2850.5 


11.93 





.4084525 


1674 


48 





.52 





.10 


0. 


19 





.03 


RRAB/EC/ESD 






203420-2508.9 


11.58 





.5262389 


666 


.44 





.90 





.17 


0. 


19 





.05 


RRAB 






203749-5735.5 


12.25 





.4199162 


-1270 


.33 





.43 





.15 


0. 


36 





.07 


RRC/EC 




* 


210741-5844.2 


13.22 





.3462376 


1479 


95 





.68 





.17 


0. 


26 





.06 


RRC/EC 




* 


211839+0612.3 


11.06 





2914601 


-1176 


.75 





.48 





.07 


0. 


15 





.01 


R.R.C 




* 


212034+1837.2 


11.50 





.5624065 


81 


.30 





.76 





.12 


0. 


16 





.03 


RRAB 






212331-3025.0 


12.29 





.3674420 


1739 


74 





.47 





.10 


0. 


21 





.03 


RRC/EC 






213826-3945.0 


12.97 





.4107031 


-1540 


.12 





.62 





.19 


0. 


30 





.09 


RRC 




* 


221556-2522.6 


11.30 





.5467383 


5 


.78 





.78 





07 


0. 


09 





.01 


RRAB 






225248-2442.2 


12.78 





.5295565 


181 


.20 


1 


.19 





.20 


0. 


17 





.05 


RRAB 






230659-4354.6 


12.66 





.2811062 


-10 


.24 





.39 





.17 


0. 


43 





.05 


RRC/EC/ESD 


BO Gru 


* 


231209-1855.4 


12.57 





.3079943 


-1349 


.89 





.47 





.13 


0. 


29 





.05 


RRC/EC 




* 


232031-1447.9 


12.46 





.6269552 


54 


.52 





.59 





.17 


0. 


29 





04 


RRAB 




* 



(Roberts et al. 1987) algorithm, with a gain of 0.5, for 200 
iterations. The frequency grid, of spacing 0.125/T, ran from 
day -1 to 10 day -1 . We then sorted the peaks and identi- 
fied pairs of peaks which were equidistant across each of the 
harmonics. Additional frequency components of the Blazhko 
periods were found surrounding higher harmonics, as well 
as a candidate for a third Blazhko period. The parame- 
ters returned by CLEAN were refined with the Levenberg- 
Marquardt algorithm according to the model 

15 

X m {t) COBpKfiit-Ti)], (1) 

i=7 

where the mean magnitude was fit to be n — 12.5980 and 
the resulting fits for the other parameters (/;, Ai, and Ti) 
are in Table 6. 



The identification and removal of these significant fre- 
quencies is illustrated in Figure 8. Panel (a) shows a Lomb- 
Scargle (Lomb 1976, Scargle 1982) periodogram of the raw 
light curve, with the first four harmonics marked with open 
triangles. A six harmonic model was fitted by least squares, 
taking observational errors into account, and subtracted 
from these data. Panels (b)-(e) show the periodogram of 
the residual light curves within 0.1 day -1 of each of the 
first three harmonics. Marked with filled triangles are fre- 
quencies corresponding to sinusoids that have already been 
subtracted, whereas marked with open triangles are fre- 
quencies identified as significant in the current panel. Pan- 
els (b), (c), and (d) correspond to the Blazhko frequen- 
cies of /si = 0.0112 day -1 , / S2 = 0.0152 day -1 , and 
f B3 — 0.0339 day -1 , matching Blazhko periods of 89.2 days, 
65.6 days, and 29.5 days, respectively. Panels (e) and (f) are 
the periodogram of the final residuals, zoomed in on the first 
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Table 3. Objects exhibiting BL2 effect. Ao is the peak-to-peak amplitude of variation as provided by ACVS and A- and A+ arc 
the peak-to-peak amplitudes of Blazhko modulation corresponding to the modulation frequency on the left and right side of the main 
pulsation frequency, respectively. See section 3.2 for details. 



ASAS ID 


V 




p 


Pbl 


A 


A_ 




A_ 


M+ 




oise 


Type 


Other ID 




[mag] 




[days] 


[days] 


[mag] 


[mag] 


[mag] 














003514-0415.0 


12.76 





.3445751 


-1616.29 





.60 





.14 


0. 


15 





96 





.03 


RRC/EC 




012848-1127.2 


12.59 





5166664 


-84.99 


1 


21 





13 


0. 


12 


1 


.06 





.01 


RR.AB /DCEP-FO 




020728-5752.2 


10.88 





3750261 


1673.36 





.44 





04 


0. 


04 


1 


.11 





.00 


RRC 


NSV00728 


021515-1048.0 


10.52 





.6234139 


112.05 





.64 





.10 


0. 


06 


1 


.55 





.01 


RRAB 


RV Cet 


030534-3116.1 


12.53 





.4964538 


-6.77 


1 


.47 





.08 


0. 


06 


1 


24 





.02 


RRAB 




031408-3446.4 


11.41 





.3124235 


-1241.77 





.52 





.06 


0. 


04 


1 


34 





.01 


RRC 




032438-2334.7 


12.09 





.6296339 


335.29 





.94 





.12 


0. 


10 


1 


26 





.01 


RRAB 




032520-6503.3 


11.24 





.4920082 


-160.64 





.88 





.11 


0. 


10 


1 


03 





.04 


RRAB / DCEP-FO: 


X Ret 


050747-3351.9 


12.05 





.4873552 


65.41 


1 


.01 





.08 


0. 


10 





77 





.03 


RRAB 


SU Col 


050747-3351.9 


12.05 





4873552 


88.98 


1 


.01 





.08 


0. 


08 





.98 





.03 


RRAB 


SU Col 


050838-5602.9 


12.10 





.5160761 


-786.91 


1 


.08 





.20 


0. 


05 


3 


79 





.03 


RRAB 


NSV 1856 


051508-4137.7 


10.45 





.4788368 


81.95 





.84 





15 


0. 


06 


2 


44 





.02 


RRAB 


RY Col 


052402-2247.4 


13.40 





6498517 


-10.23 


1 


19 





10 


0. 


11 





.94 





.03 


RRAB 




052840-5316.2 


13.38 





.3678062 


-453.43 





.67 





.20 


0. 


17 


1 


.21 





.05 


RRC 




053628-3837.0 


12.64 





3714727 


-1394.31 





.57 





.09 


0. 


14 





.61 





.01 


RRC 




054843-1627.0 


12.88 





.3767273 


1663.06 





.88 





.15 


0. 


17 





93 





.03 


RRAB/RRC: 




070001-3732.5 


11.85 





.4941843 


116.96 





.89 





.13 


0. 


11 


1 


21 





.02 


RRAB 


NSV03331 


080249-5913.5 


11.75 





3541891 


-1185.26 





.28 





.06 


0. 


06 





90 





.01 


RRC/EC/ESD 


_ 


080318-2530.1 


10.86 





.2697723 


-1542.02 





.43 





.04 


0. 


03 


1 


30 





.01 


RRC 


NSV03882 


085448-8317.0 


12.27 





.4778598 


-36.79 


1 


.30 





.14 


0. 


10 


1 


.40 





.03 


RRAB 


NSV04350 


091349-0919.1 


10.71 





.5372290 


26.30 





.90 





.15 


0. 


16 





97 





.03 


RRAB /DCEP-FO 


SZ Hya 


094438-4552.6 


11.58 





5735076 


66.35 





.95 





.05 


0. 


06 


(1 


.98 





.01 


RRAB 


CD Vel 


103203-3010.6 


11.56 





3304459 


1731.60 





.47 





.06 


0. 


06 


1 


.00 





.01 


RRC 


NSV04885 


103246-3423. 1 


13.39 





.2448541 


-277.05 





.78 





.12 


0. 


18 





65 





.04 


RRC 




105303-4954.4 


10.76 





5274139 


58.68 





85 





.08 


0. 


07 


1 


17 





.02 


RRAB 


AF Vel 


110522-2641.0 


11.68 





.2944559 


-7.40 





.38 





.08 


0. 


07 


1 


.14 





.02 


RRC 




112027-4338.8 


11.16 





.3795948 


1480.38 





.28 





.03 


0. 


04 





92 





.03 


RRC/EC/ESD 




114555-5922.7 


11.16 





.4531949 


-79.01 


1 


19 





.07 


0. 


08 





.90 





.02 


RRAB 


BI Ccn 


120447-2740.7 


9.76 





.6503243 


71.79 





.62 





.13 


0. 


08 


1 


.61 





.03 


RRAB 


IK Hya 


120942-3457.4 


13.20 





3436271 


1611.60 





.66 





.15 


0. 


10 


1 


.47 





.03 


RRC 


EL Hya 


121206-2612.8 


12.65 





3987747 


48.46 


1 


31 





.20 


0. 


15 


1 


.27 





.04 


RRAB 




123030-2602.9 


9.81 





.4785475 


63.29 





.94 





05 


0. 


03 


1 


.41 





.01 


RRAB 


SV Hya 


132922-0553.0 


11.24 





.5763497 


1398.99 





.87 





.14 


0. 


08 


1 


76 





.04 


RRAB 




140324-3624.4 


10.71 





.4939666 


-1650.71 





.98 





.15 


0. 


17 





86 





.04 


RRAB 


V0674 Cen 


141345-2254.7 


12.05 





.4479434 


26.50 


1 


.00 





.30 


0. 


18 


1 


.65 





.04 


RRAB 




150924-4319.6 


12.47 





.3821508 


42.49 


1 


.15 





.26 


0. 


12 


2 


08 





.04 


RRAB 


FU Lup 


151849-1000.0 


12.03 





.3364272 


802.95 





.50 





.17 


0. 


03 


4 


84 





.08 


RRC 




153517-2420.2 


11.24 





3067770 


1560.06 





.54 





.06 


0. 


06 


1 


.10 





.02 


RRC 


CG Lib 


153830-6906.4 


12.26 





.6224747 


-118.05 





.89 





11 


0. 


08 


1 


36 





.02 


RRAB 




155552-2148.6 


11.38 





.2541338 


1699.52 





.46 





10 


0. 


09 


1 


.12 





.02 


RRC 




160204+1728.8 


10.72 





2308114 


11.49 





.44 





04 


0. 


05 





.82 





.04 


RRC 


LS Her 


160204+1728.8 


10.72 





.2308114 


12.76 





.44 





07 


0. 


10 





74 





.04 


RRC 


LS Her 


162158+0244.5 


12.47 





3238044 


-8.11 





.52 





12 


0. 


12 





.99 





03 


RRC/EC 




162811+0304.3 


13.09 





.5970104 


-26.28 


1 


.65 





.14 


0. 


15 





92 





.03 


RRAB 




163225-8354.2 


8.99 





5425804 


-143.73 





.85 





.08 


0. 


08 





.99 





.04 


RRAB 


UV Oct 


172932-5548.3 


12.89 





.3273052 


-1610.05 





.68 





.20 


0. 


17 


1 


.19 





.02 


RRC 


EZ Ara 


174048-3132.6 


10.60 





4272967 


-504.03 





.85 





.05 


0. 


04 


1 


.17 





.02 


RRAB 


V0494 Sco 


174202-4633.7 


10.72 





3115788 


-1706.78 





.49 





.10 


0. 


10 


1 


.01 





.03 


RRC 




175911-4926.0 


10.11 





.4518593 


-49.37 





.94 





.08 


0. 


08 





97 





.02 


RRAB 


S Ara 


183441-6527.0 


11.88 





.4769536 


-173.70 


1 


17 





11 


0. 


09 


1 


.22 





.02 


RRAB 


BH Pav 


184758-3744.4 


10.26 





.5893445 


-59.96 





.59 





.02 


0. 


02 


1 


36 





.01 


RRAB 


V0413 CrA 


192824-1852.4 


12.57 





.3563567 


1572.33 





.61 





17 


0. 


19 





93 





.05 


RRC 


" 


193538-7409.9 


12.51 





.3499993 


1608.49 





.60 





11 


0. 


12 





.97 





.02 


RRC/EC 




195142-6244.1 


11.86 





.5514395 


557.17 


1 


.14 





.12 


0. 


11 


1 


10 





.02 


RRAB 


FO Pav 


195927-3400.1 


11.88 





.3797209 


45.69 





.77 





.14 


0. 


11 


1 


.26 





.04 


RRAB 




200910-4149.5 


12.61 





.4401943 


45.39 


1 


15 





09 


0. 


08 


1 


.16 





.02 


RRAB 


V2239 Sgr 


202044-4107.1 


10.93 





5529452 


1331.74 





.86 





12 


0. 


10 


1 


.27 





.02 


RRAB 


V1645 Sgr 


203145-2158.7 


11.25 





3107152 


792.83 





.39 





.04 


0. 


04 


1 


.11 





.01 


RRC/EC 




204440-2402.7 


12.75 





.2053330 


-6.64 





.36 





.06 


0. 


04 


1 


26 





.01 


RRC/DSCT/EC 




210129-1513.8 


10.47 





.4477465 


231.66 


1 


.04 





.17 


0. 


07 


2 


35 





.02 


RRAB 


RV Cap 


212433-5712.1 


12.95 





.6051401 


-133.38 


1 


.14 





.14 


0. 


09 


1 


r>4 





.02 


RRAB 




214101+0109.6 


12.43 





.6156709 


-522.58 





.57 





.13 


0. 


10 


1 


.36 





.02 


RRAB 




214719-8739.1 


12.34 





.4580005 


244.20 


1 


.25 





.09 


0. 


09 


1 


.01 





.04 


RRAB 


RS Oct 


215336-8246.8 


11.34 





.6218493 


144.12 





.98 





.04 


0. 


05 





.80 





.01 


RRAB 


SS Oct 


220254-2131.5 


10.73 





.3637073 


1413.43 





.46 





07 


0. 


07 


1 


04 





.01 


RRC 


BV Aqr 


222539-0756.5 


10.66 





.4052637 


1618.65 





.37 





07 


0. 


07 


1 


.10 





.01 


RRC 


GP Aqr 


225131-3006.2 


13.10 





.3384769 


-1681.80 





.62 





12 


0. 


12 





.99 





.03 


RRC/EC 




225323+0846.1 


12.60 





.4930493 


-348.58 


1 


.51 





18 


0. 


21 





.88 





.04 


RRAB 




225518-2317.6 


12.99 





.3935794 


1557.88 





.62 





13 


0. 


10 


1 


.29 





.01 


RRC/EC 




233951-1644.4 


12.16 





.3553741 


875.96 





.52 





17 


0. 


13 


1 


26 





.06 


RRC 
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Table 4. Objects exhibiting long period changes (PC). See section 3.2 for details. 





V 
[mag] 




p 
[days] 


[days] 


[mag] 


A 1 

[mag] 




/J. 




oise 


Type 


Other ID 


010117-4556.6 


11 


70 





3738269 


2085.07 


0.27 


0.08 





29 





.02 


RRC 




011831-6755.1 


11.42 





.4057948 


-1702.13 


0.49 


0.10 





.21 





.04 


RRC 


AM Tuc 


032337+1540.0 


9.66 





.2384770 


-1738.83 


0.13 


0.02 





17 





.00 


RRC/EC/ESD 




033108+0713.4 


10 


.64 





5281963 


1116.69 


0.18 


0.04 





.24 





.02 


RRAB/EC/ESD 




040500-4457.1 


12 


.66 





.5668385 


-1533.98 


0.78 


0.05 





.06 





.01 


RRAB 




045213-5620.8 


13 


.59 





.3503107 


2050.02 


0.75 


0.21 





28 





.06 


RRC 




052406-6925.2 


13 


.18 





.5531259 


2668.80 


2.21 


0.18 





08 





.02 


RRAB 


TY Dor 


052840-5316.2 


13 


.38 





.3678062 


2178.17 


0.67 


0.16 





.24 





.05 


RRC 




062838-3848.5 




.53 





.4834677 


2282.06 


1.41 


0.10 





.07 





.03 


RRAB 




065307-5935.7 


11 


.47 





.7370705 


-2123.14 


0.85 


0.05 





.06 





.02 


RRAB 


IU Car 


085816-3022.1 


10 


.60 





5107257 


-2403.27 


0.12 


0.03 





.22 





.01 


RRAB/DCEP-PO/EC/ESD 




104749-0308.8 


13 


.04 





.3632554 


2236.14 


1.02 


0.14 





13 





.06 


R.R.AB/DCEP-FO/EC 




112027-4338.8 


11 


.16 





.3795948 


-2346.87 


0.28 


0.04 





14 





.03 


RRC/EC/ESD 




112715-2510.4 


13 


.39 





.5768075 


-2457.61 


1.66 


0.15 





.09 





.06 


RRAB/DCEP-FO 


- 


113701-0600.4 


12 


.57 





.4022837 


1880.76 


0.60 


0.09 





.15 





.01 


RRC/EC 




140324-3624.4 


10 


.71 





.4939666 


-1650.71 


0.98 


0.17 





18 





.04 


RRAB 


V0674 Ccn 


151849-1000.0 


12 


.03 





.3364272 


-2192.98 


0.50 


0.03 





07 





.08 


RRC 




151951-0950.0 


13 


.38 





.3235636 


-8591.07 


0.70 


0.00 





.00 





.02 


RRC: 




154504+1736.7 


12 


.82 





.5736250 


1309.24 


1.27 


0.09 





.07 





.07 


RRAB 




161301-2406.9 


12 


.41 





.6024391 


2036.66 


0.43 


0.07 





17 





.02 


RRAB: 




195927-3400.1 


11 


.88 





.3797209 


-2285.19 


0.77 


0.12 





.16 





.04 


RRAB 




210741-5844.2 


13 


.22 





.3462376 


-1839.25 


0.68 


0.17 





.25 





.06 


RRC/EC 




212045-5649.2 


12 


.38 





.3680388 


2286.76 


0.26 


0.05 





.20 





.02 


RRC/EC 




212756-6702.6 


12 


.73 





.3449728 


2410.80 


0.43 


0.09 





.20 





.03 


RRC 




220454-6635.0 


12 


.05 





.5640177 


1722.65 


0.88 


0.14 





.15 





.05 


RRAB 


NSV14009 


221039-5049.8 


12 


.32 





.3306480 


-1720.28 


0.70 


0.15 





.21 





.05 


RRC 




224131-0628.6 


11 


.86 





.5744469 


-1771.17 


0.73 


0.15 





.21 





.02 


RRAB 


HH Aqr 


225559-2709.9 


12 


.71 





.3103066 


-2197.80 


0.37 


0.08 





.22 





.02 


RRC/EC 




232342-4157.4 


12 


.72 





.4424689 


-2144.08 


0.33 


0.07 





21 





.01 


RRC/EC/ESD 





Table 5. RR Lyrac Blazhko stars with two pairs of BL2s. See section 3.2 for details. 



ASAS ID 


V 


^0 


P 


Pbl 




A 2 


A_/A + 


Type 


Other ID 




[mag] 


[mag] 


[days] 


[days] 


[mag] 


[mag] 








050747-3351.9 


12.05 


1.01 


0.4873552 


65.8 


0.08 


0.10 


0.77 


RRAB 


SU Col 










89.3 


0.08 


0.08 


0.98 






160204+1728.8 


10.72 


0.44 


0.2308114 


11.5 


0.04 


0.05 


0.82 


RRC 


LS Her 










12.8 


0.07 


0.10 


0.74 







three harmonics and of the entire interval, respectively. The 
two sets of BL2 frequencies near the main peak — left hand 
side of panel (b) — is what led us to investigate this star in 
detail. In some cases apparently significant peaks disappear 
upon the removal of another peak; this is a symptom of alias- 
ing. An example can be clearly seen when the peak at 4.04 
day -1 in panel (c) disappears upon removal of the peak at 
/n = 2.037 day -1 . We used the CLEAN algorithm to avoid 
precisely this type of confusion. 

Despite our efforts using CLEAN, the reality of the 
third Blazhko period is not conclusively established because 
there is a peak in the window function at /' = 2.0027 day -1 
which couples fs2 to /b3- The following numerical relation- 
ships approximately hold: 

/l5-/' « fl2 (2) 
/l4 + /' W /2-/S2. (3) 

Equation 2 implies /15 could be a nearly 2 day -1 alias of 
/12. In addition, we note that the fitted parameters of these 
sinusoids have a large covariance, about half the geometric 
mean of the individual parameters' variances. Equation 3 
gives the reason for this covariance, although our analysis 
did not identify a significant sinusoid at fa — /s2; if it exists, 
it is apparently buried in noise. 



The status of SU Col as an RR Lyrae variable was es- 
tablished by Gessner (1985). Berdnikov and Turner (1996) 
presented an O-C diagram and refined the ephemeris of max- 
imum light, giving: 

Max HJD = 2450363.870(2) + 0.48735842(11) x E. (4) 

The uncertainties are in parentheses, indicating the pre- 
cision of the final and final two digits for epoch and period, 
respectively. 

The maximum with number E — 5008 is near the 
center of our dataset, and gives a time of maximum C = 
2452804.561. Although the cadence of our dataset does not 
allow individual maxima to be found, the nearby maximum 
of the six harmonic model is at O = 2452804.558 ± 0.010, 
matching the published ephemeris to within the error (O — 
C = 0.002 ±0.010). Our best-fitting main period of 0.487355 
d ±2 x 10~ 6 d, also marginally fits the ephemeris of equation 
4. 

The Hipparcos and Tycho Catalog (Perryman and ESA 
1997) reported SU Col as an RR Lyrae as well, reporting 
the epoch of maximum as Oh = 2448500.2120 ± 0.0010. 
The uncertainty was only reported to order-of-magnitude, 
and is based on a three harmonic model; judging by the 
folded light curve it is probably at least a factor of a few 
too low. This corresponds to maximum number E = —3825, 
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Table 6. Best-fitted parameters for the detected frequencies in 
the Blazhko star SU Col. The parameters listed below refer to 
equation 1. Times of light maximum are listed as HJD-2450000. 
Sec section 3.3. 



/ (day- 1 ) A (mag) T (day) note 



2.0518901 


0.69 


2804.6078 






h 




4.1037803 


0.30 


2804.5774 




h 




2/i 


6.1556702 


0.20 


2804.5569 




h 




3/i 


8.2075605 


0.11 


2804.5484 




h 




4/i 


10.259451 


0.08 


2804.6462 




h 




5/i 


12.311340 


0.05 


2804.6246 




h 




6/i 


2.0407786 


0.08 


2804.6775 


h 




h 


— /si 


2.0631858 


0.09 


2804.4963 


h 




h 


+ /si 


4.0926081 


0.06 


2804.6302 


h 




h 


— /si 


4.1149220 


0.06 


2804.5249 


ho 




h 


+ /si 


2.0366963 


0.09 


2804.5442 


hi 




h 


- }B2 


2.0671779 


0.09 


2804.5708 


hi 




h 


+ /S2 


6.1403548 


0.06 


2804.5370 


fl3 




h 


— /S2 


2.0857132 


0.05 


2804.6723 






h 


+ /S3 


4.0697479 


0.06 


2804.5647 


fit 




h 


~ /S3 



whose calculated value is Ch = 2448500.2114, giving O — 
C = 0.0006 ± 0.0022. Therefore equation 4 satisfies three 
modern data sets over a span of about fifteen years, despite 
the complication that the Blazhko effect causes on the time 
of maximum light. 

The status of SU Col as a Blazhko star has not been 
previously recognized. Wils, P. and Sodor (2005) analyzed 
the ASAS RR Lyrae stars for Blazhko behavior, although 
we suspect that SU Col was not analyzed because its folded 
light curve looks like it is affected by noise rather than 
the Blazhko effect, and visually inspecting the folded light 
curves was their first cut. Kovacs (2005) also identified 
Blazhko stars in ASAS, however, he restricted his attention 
to stars whose metallicity was already measured. 

Recently it was discovered that another star with multi- 
ple Blazhko periods, namely XZ Cyg, is changing its Blazhko 
period over time (LaCluyze et al. 2004). The authors dis- 
cover that these changes are anticorrelated with observed 
changes in the primary period of this star. We looked for 
a possible Blazhko period change for SU Col. In particular, 
we divided the light curve in two, but each section had the 
same period. Also, an O-C diagram for this object does not 
show any long term primary period changes (although per- 
turbations by the Blazhko effect are seen), so the multiple 
Blazhko periods of SU Col cannot be explained analogously 
to XZ Cyg behaviour. 

Future observations with a high cadence, spanning a few 
hundred days, would be extremely valuable for understand- 
ing the multiply periodic nature of the Blazhko effect in this 
star. Such a campaign has recently reported similar double 
Blazhko behavoir for UZ UMa (Sodor et al. 2006). How- 
ever, a campaign with telescopes distributed in longitude 
may be necessary to resolve all aliases. SU Col is consider- 
ably brighter than the class of objects with two sets of BL2 
peaks found by Collinge et al. (2006), so even time-resolved 
spectroscopy is not prohibitive. Discovering rare objects that 
are bright enough to follow up is one of the strengths of wide- 
field surveys like ASAS, and we urge follow-up observations 
of this object. 



SU Col (ASAS 050747-3351.9) 
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Figure 8. An iterative subtraction of extra frequency compo- 
nents of SU Col illustrates the rich multiperiodic behaviour of 
this Blazhko star, (a) The Lomb-Scarglc periodogram of the raw 
light curve, (b) Periodogram after removing frequency sinusoids 
at the main frequency and 5 higher harmonics. The three panels 
are zoomed in on the region surrounding the first three harmonics. 
A number of Blazhko peaks are noticeable, (c, d, e) Periodograms 
after removing the extra peaks associated with /si, /s2i and 
/b3 respectively (sec tabic 7). (f) The final periodogram after 
subtracting the 15 peaks listed in Tabic 6. In each panel, open 
triangles identify the peaks that are most significant, and filled 
triangles identify frequencies at which peaks have already been 
subtracted. See section 3.3. 

4 DOUBLE-MODE RR LYRAE 

RRd variables are objects pulsating in two radial modes: 
the fundamental mode and the first overtone. Characteris- 
tic is the period ratio of these modes, Pi/Pq ~ 0.744, where 
Po and Pi correspond to the periods of fundamental mode 
and the first overtone respectively. So far, 24 such isolated 
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objects have been detected in the field of our Galaxy, 17 
of them are listed by Wils (2006, see references therein) 
and 7 by Bernhard & Wils (2006), without including fainter 
objects of the Galactic Bulge (Mizerski, 2003; Moskalik & 
Poretti, 2003; Pigulski et al., 2003), and the foreground stars 
of the dwarf galaxy in Sagittarius (Cseresnjes, 2001); 13 of 
them were found in the ASAS data. 

Very useful in the study of RRd stars is a Petersen 
diagram (Petersen 1973), on which the fundamental mode 
period (Po) is plotted against the ratio of the overtone and 
fundamental mode periods (Pi/Pq). Several attempts were 
made to constrain theoretically the region occupied by RRd 
stars on a Petersen diagram (Bono et al. 1996, Popiclski 
et al. 2000, Bragaglia et al. 2001, Kovacs 2001). Most re- 
cently Szabo et al. (2004) combined both the double mode 
pulsational and evolutional theories to construct such a the- 
oretically permitted region. We searched a slightly bigger 
region (0.40 < P < 0.65 and 0.735 < Pi/P„ < 0.755) for 
RRd variables without putting any additional constraints. 
The condition that a star had to fulfill to classify as a RRd 
candidate was that the ratio of the two highest amplitude 
frequencies falls within the desired period ratio range. This 
procedure provided us with 35 candidates. 

We visually examined those and 16 turned out to be 
spurious detections. So in our search for multiperiodicity in 
ASAS data we detect 19 double mode RR Lyrae variables, 
among which 4 are newly identified RRd stars. This sums 
up to 17 RRd stars altogether discovered in ACVS, out of 28 
objects now known in the field of our Galaxy (on both south- 
ern and northern hemisphere) , giving the ASAS project the 
majority of discoveries. The light curves of newly discov- 
ered objects, phased with both pulsation periods are given 
in Figure 9. 

Surprisingly, we did not detect 4 objects that were dis- 
covered as RRd before and should be present in the ASAS 
database. After a closer investigation of these objects we 
found that two of them, namely ASAS 030528-3058.7 (Bern- 
hard and Wils, 2006) and CU Com = ASAS 122447+2224.5 
(Clementini et al., 2000) are phased in ACVS with one of 
the aliases of the true period, so the period ratio fell outside 
the searched region. For the third one (Jerzykiewicz and 
Wenzel, 1977), namely AQ Leo = ASAS 112355+1019.0, 
we did not find a significant peak among 50 highest that 
would correspond to the second pulsation period. However, 
the light curve phased with this period has a modulation 
that appears significant, which apparently suggests that the 
CLEAN algorithm failed in this case. The fourth case was 
an object found by Wils et al. (2006) in NSVS data in an 
overlap region, namely V458 Her = ASAS 170831+1831.3. 
We do not confirm any significant secondary frequency in 
this light curve with the longer dataset. 

The incidence ratio, which we define as a number of 
double mode pulsators divided by the number of RRc stars 
is RRd/RRc = 0.025. This value is small, comparing to the 
one determined by Alcock et al. (2000) for the LMC, which 
is 0.134. The ratio for the Galactic bulge calculated by Miz- 
erski (2003) is 0.007. 

Parameters of these objects, such as their coordinates, 
magnitudes, pulsation periods and amplitudes, are listed in 
Table 7. Also another star ID is given if applicable and a 
reference to the discovery paper. For stars that are new dis- 
coveries, there is a "*" instead of a reference. One of the 



ASAS 032820-6458.7 PO - 0.477 d 



ASAS 032820-6458.7 P1 = 0.355 d 



13.1 
13.2 



!."'■ ■ V,ri".*i-s ;* :--J.;S\tS-y:-. ■- 
^'■.VS-«?H;>.Mf^'-.VS-«?H;>.y 



■:m... 



ASAS 040054-4923.8 PO . 0.559 d 



14.4 
14.6 



ASAS 040054-4923.8 P1 . 0.41 7 d 



0.5 1 



ASAS 1 33439+241 6.6 P0 = 0.488 d 



ASAS 133439+2416.6 P1 -0.363 d 



1 12 



ASAS 211848-3430.4 PO . 0.505 d 



ASAS 21 1848-3430.4 P1 . 0.376 d 



13.5 
13.6 
13.7 



•<-j>./: • >-i >..".. 



Figure 9. Four newly discovered double-mode pulsators are rep- 
resented in 4 rows. Each row contains a light curve phased with 
the main pulsation period (left) and the first overtone period 
(right). See section 4 for details. 



newly found double mode stars, namely 032820-6458.7 (or 
SW Ret), is a secure detection based on its period ratio 
of 0.744417, but it is classified in GCVS as W UMa type 
eclipsing variable. However, due to slope asymmetry, ACVS 
unambiguously classifies this object as an RRc star. 

It is known that the majority of RRd stars have a first 
pulsation mode with a lower amplitude than the first over- 
tone (eg. Alcock et al. 2000). The ASAS detections confirm 
this: 18 out of 19 objects have amplitude ratio Ai/Ao > 1. 

Figure 10 shows a Petersen diagram for double-mode 
RR Lyrae in ASAS data. All stars follow the sequence ob- 
served up to date by many RR Lyrae (see Figure 13 in 
Clementini et al. 2004). The new discoveries on our diagram 
are plotted with solid circles and known objects with open 
circles. 



5 SUMMARY 

We performed a multiperiodicity search in the ASAS Galac- 
tic field RR Lyrae data, consisting of 1435 RRab and 756 
RRc stars. It resulted in finding 160 objects, 122 of them ex- 
hibiting Blazhko effect (BL), 29 long period changes (PC), 
and 19 double-mode behaviour (RRd stars). 

There are 73 BL stars inside the RRab group, and 49 in 
RRc, giving a total of 5.6% BL in the whole sample, which 
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Table 7. Parameters for double-mode RR Lyrae stars (RRd) found in ASAS data. See section 4. 



ASAS ID 


Y 


1 


An 


p. 
1 1 


A i 


Pi / Pn 


At / An 


Other ID 




(RA DEC) 


[mag] 


[uay sj 


[mag] 


T rl a ^7 c 1 
[oay sj 


[mag] 










032820-6458.7 


12.56 


0.4766242 


0.07 


0.3548111 


0.18 


0.7444253 


2.73 


SW Ret 


* 


040054-4923.8 


13.32 


0.5585883 


0.10 


0.4166705 


0.16 


0.7459348 


1.61 




* 


081610-6644.8 


12.36 


0.5172180 


0.14 


0.3852143 


0.21 


0.7447813 


1.47 


GSC 8936-2145 


(WIls & Otero 2005) 


084747-0339.1 


10.53 


0.5639063 


0.08 


0.4207890 


0.20 


0.7462037 


2.59 


GSC 4868-0831 


(WIls & Otero 2005) 


122509-2139.9 


12.19 


0.5445154 


0.03 


0.4064419 


0.20 


0.7464286 


6.01 




(Bernhard & Wils 2006) 


133439+2416.6 


12.44 


0.4879017 


0.12 


0.3630583 


0.23 


0.7441218 


1.85 


BS Com 


(Bragaglia ct al. 2003) 


141539+0010.1 


13.35 


0.4819318 


0.12 


0.3588418 


0.18 


0.7445904 


1.57 




* 


151735-0105.3 


10.96 


0.4713494 


0.15 


0.3507304 


0.21 


0.7440985 


1.41 


V0372 Scr 


(Garcia-Mclendo et al. 2001) 


173726+1122.4 


12.76 


0.4633236 


0.17 


0.3442391 


0.17 


0.7429777 


1.02 


V2493 Oph 


(Garcia-Melendo & Clement 1997) 


183952-3200.9 


11.81 


0.4612583 


0.20 


0.3424752 


0.18 


0.7424802 


0.89 


GSC 7411-1269 


(WIls & Otero 2005) 


184035-5350.7 


12.21 


0.4790651 


0.07 


0.3563601 


0.10 


0.7438656 


1.49 




(Bernhard & Wils 2006) 


193933-6528.9 


12.89 


0.4915148 


0.17 


0.3657414 


0.20 


0.7441106 


1.20 


GSC 9092-1397 


(Wils 2006) 


195612-5043.7 


12.12 


0.4678191 


0.13 


0.3477846 


0.15 


0.7434169 


1.17 


GSC 8403-0647 


(WIls & Otero 2005) 


210726+0110.3 


13.06 


0.4772373 


0.17 


0.3549859 


0.22 


0.7438352 


1.27 




(Bernhard & Wils 2006) 


211848-3430.4 


12.84 


0.5048597 


0.10 


0.3763662 


0.22 


0.7454867 


2.12 




* 


212721-1908.0 


13.06 


0.4730220 


0.13 


0.3520623 


0.19 


0.7442830 


1.46 


NSV13710 


(Bernhard & Wils 2006) 


213437-4907.5 


12.00 


0.4879954 


0.14 


0.3631872 


0.19 


0.7442432 


1.38 


Z Gru 


(Wils 2006) 


230449-3345.3 


12.99 


0.4978410 


0.14 


0.3710152 


0.21 


0.7452483 


1.47 




(Bernhard & Wils 2006) 


235622-5329.4 


12.58 


0.5668088 


0.10 


0.4224954 


0.20 


0.7453932 


2.05 


NSV14764 


(Bernhard & Wils 2006) 
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Figure 10. A Petersen diagram for galactic field RRd variables 
discover during multiperiodicity search on ASAS data. Known 
objects arc plotted with open circles, while newly identified stars 
are drawn with filled circles. See section 4 for details. 

is few times less than the percentage observed in Magel- 
lanic Clouds and Galactic Bulge, but similar to the number 
found in previous searches in Galactic field survey catalogs. 
The analysis of the Blazhko sample reveals a significant dif- 
ference between the behaviour of RRab and RRc groups. 
While RRab stars have a wide, unimodal distribution in the 
Blazhko period, RRc stars occupy either the short Blazhko 
period region around 10 days, or a long Pbl area, starting 
around 300 days, with a condensation roughly around 1500 
days. There is a lack of Blazhko stars among high period 
(P > 0.65 d) RRab stars, which may give insight into the 
physics of the Blazhko effect, but we suspect simply indi- 
cates the misclassification of RRab stars. We do not con- 
firm any correlation between the pulsation period and the 
Blazhko period. The absence of visible differences between 
the BL1 and BL2 groups suggests that the same mechanism 
is responsible for the Blazhko effect in both cases. 



During our study we discovered the Blazhko effect with 
multiple periods (BL2x2) in object ASAS 050747-3351.9 = 
SU Col. were identified with periodogram peaks near the 
first three harmonics of the main pulsation. 

We also identified 19 double-mode RR Lyrae, 4 of them 
being new discoveries, one of these (SW Ret) being reclassi- 
fied from an eclipsing contact binary to a RR Lyrae variable. 
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